Discrete element model for general polyhedra

Abstract We present a version of the Discrete Element Method considering particles as rigid polyhedra. The Principle Virtual Work is employed basis for multibody dynamics model. Each particle surface split into sub-regions, which are tracked contact with other sub-regions neighboring particles. Contact interactions modeled pointwise, vertex-face, edge-edge, vertex-edge and vertex-vertex interactions. General polyhedra triangular faces considered particles, permitting multiple pointwise automatically detected along model evolution. propose combined interface law composed penalty barrier approach, to fulfill constraints. Numerical examples demonstrate that can handle normal frictional effects in robust manner. These include simulations convex non-convex showing potential applicability materials complex shaped such sand railway ballast.